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Abstract 

We argue that the coarse-grained dynamics of the zero-range process in the condensation regime 
can be described by an extension of the standard hydrodynamic equation obtained from Eulerian 
scahng even though the system is not locahy stationary. Our result is supported by Monte Carlo 
simulations. 
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I. INTRODUCTION 

The zero-range process (ZRP) was introduced in 1970 by Spitzer system of 

interacting random walks, where each lattice site k is occupied by ra^. particles which hop 
randomly to other sites. The hopping rates Wn depend only on the number of particles n 
at the departure site. Under certain conditions on the rates Wn and the particle density 
(see below) the grand-canonical stationary distribution is a product measure, i.e. there are 
no correlations between different sites [2]. An exact large-scale description of the dynamics 
has been proved for arbitrary initial densities in terms of a hydrodynamic equation for the 
coarse grained particle density p{x,t'), provided the rates are non- decreasing, i.e., Wn+i > 

m 

. In this so-called attractive case the density satifies the continuity equation 

dt'P + d,j{p) = (1) 

where j{p) is the stationary current in the grand-canonical distribution with density p. 

Depending on the choice of rates, in non-attractive systems a rich and rather varied 
dynamical and stationary behaviour emerges, for a recent review see jsl]. In particular, the 
model may admit a condensation phenomenon analogous to Bose-Einstein condensation. 
In a periodic chain one then finds that, above a critical density pc, a finite fraction of all 
particles in the system accumulate at a randomly selected site, whereas all other sites have 
an average density pc js, 3, S Q]- T 
in terms of a coarsening process 



al. 



le late-time dynamics of condensation has been studied 
[lo| . For analysis of the early time range Kaupuzs et 



ll| have generalized an equivalent exclusion process for traffic flow They observed 
a metastable regime as a precursor of the coarsening process. The metastable state with a 
current above the critical current persists for some finite time range which depends on the 
hopping rates. 

If valid, the hydrodynamic description ([1]) would apply to times much later than this 
metastable regime. Interestingly, however, the proofs for the hydrodynamic limit do not work 
when the condition on the rates Wn for condensation is met. The reason for this failure is not 
a minor technical issue of the standard hydrodynamic approach but lack of local stationarity 
during the coarsening process, the early stages of which fall within the hydrodynamic time 
regime. This profound violation of one of the most basic assumptions of hydrodynamic 
theory in conjunction with the highly discontinuous condensation phenomenon may lead 
one to suspect that ([1]) might not be valid for initial density profiles where some region of 



space is above the critical density pc- It is the aim of this paper to show that such a view, 
even though well- motivated, is overly pessimistic. We argue that, properly interpreted, the 
hydrodynamic limit ([T]) is robust and valid also in the condensation regime. The theoretical 
analysis of Sec. Illll is supported by Monte Carlo simulation in Sec. IIVI 

II. TOTALLY ASYMMETRIC ZERO-RANGE PROCESS 

To set the stage for the new ideas in the next section, we now define the details of 
our model and review some known results including the hydrodynamic limit for subcritical 
densities. 

For definiteness we consider here the one-dimensional totally asymmetric zero-range pro- 
cess (TAZRP) with periodic boundary conditions. The one-dimensional ZRP has attracted 
particular interest since it can be mapped to an exclusion processes for which double- 
occupancy of sites is forbidden. The ZRP-particles are turned into particle clusters between 
consecutive vacant sites (which correspond to the sites on which the ZRP is defined). Con- 
densation thus corresponds to phase separation between a macroscopic particle cluster and a 
disordered domain which also contains vacant sites. The n-dependence of the hopping rates 
corresponds to a length-dependent rate of detachment of a particle from a cluster of length 
n. The ZRP has thus served for deriving a quantitative criterion for the existence of non- 
equilibrium phase separation 13| in the otherwise not yet well-understood driven diffusive 

nn 

systems with two conservation laws |14J. llSj. Alternatively one can map ZRP-particles onto 
strings of vacant sites between consecutive particles (which correspond to the sites on which 
the ZRP is defined). In this mapping the n-dependence of the hopping rates Wn translates 
into a distance-dependent hopping rate for exclusion particles as may be expected for one- 
dimensional driven motion which is embedded in three-dimensional space. This could be of 
importance not only for the understanding of vehicular traffic, but also for obtaining insight 
into the dynamics of ribosomes along m-RNA [16] or the motion of molecular motors along 
microtubuli 17]. Other fields of application of the one- dimensional ZRP include experiments 

nn 

on condensation and metastability in granular media ^8|, |19| . 



If P ^ Pc then the grand-canonical stationary product measure has one-site marginals 
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P*{n) = Prob[nfc = n] given by 

nn) = ^rflw-\ (2) 

i=l 

Here the empty product (n = 0) is defined to be 1, 

oo n 

z = J2rll^;' (3) 

n=0 4=1 

is the local "partition function", and is the fugacity which determines the density p = 
(p {d /d(j)) In Z {(/)). Due to particle conservation the product distribution defined by ([2]) is 
stationary for every for which Z exists. In the TAZRP particles hop with rate Wn from 
site k to site /c + 1 on a periodic chain with L sites. This process satisfies pairwise balance 
{2^, leading to a macroscopic stationary current 

00 

j = J2wnP*{n) = <p. (4) 

n=X 

The convexity of Z ensures that the current is an increasing function of the density. In the 
case of condensation the radius (pc of convergence of the partition function is finite, with a 
critical density pc < 00 as approaches 0c- The product measure does not exist for densities 
beyond pc- 

An intuitively convenient starting point for a coarse grained hydrodynamic description 
of the dynamics of the TAZRP is the lattice continuity equation 

j^Pk{t) = jk-i{t) - 3k{t) (5) 

for the expected local density pkif) = {nk{t)), starting from some initial distribution. Here 

00 

jk{t) = Wn,P{nk,t) (6) 

is the expected local current with the probability P{nk, t) = { Sn,.{t),nk ) of finding particles 
on site k at time t. We consider Eulerian scaling where the lattice constant a (so far implicitly 
assumed to be unity) is taken to zero and the system is studied for rescaled time t' = ta, i.e., 
the microscopic time t is taken to infinity such that the macroscopic time t' is fixed. Since we 
are working with a periodic chain with L sites we take a = 1/L and correspondingly t = Lt'. 
In the hydrodynamic limit L — 00 the discrete chain of L sites becomes a continuous ring 
of circumference 1. 
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We first consider the subcritical regime where initially pk < Pc everywhere on the lattice. 
In this case one obtains formally dt'p{x,t') + dxj{x,t') = by setting k = xL in the lattice 
continuity equation and Taylor expanding in 1/L. In order to arrive at the continuity 
equation ([T]) one proves local stationarity which can be done rigorously for the ZRP and 
other lattice gas models under fairly generic circumstances. Local stationarity means that 
in the local environment of the point x, i.e., in a large but finite lattice region around the 
lattice point x = kL, the system is found in its stationary state. Physically, it follows from 
considering the limit where the microscopic time t oo, which allows all nonconserved (fast) 
local degrees of freedom to relax to their local stationarity distribution at the local density 
p (which because of the conservation law is a slow dynamical variable). The identification of 
the expected local density p{x,t') (appearing in (jS])) with the coarse grained density of the 
ZRP (appearing in ([1])) comes from the law of large numbers and local stationarity ensures 
that j{x,t') is the stationary current j{p{x,t')) computed from the product measure; for 



details see 



2l|. 



The hydrodynamic equation ([T]) can be solved by writing dfp + dpj{p)dxP = and using 
the method of characteristics. These are the lines x{t') = Vchait' along which the density 
remains constant. The characteristic velocity is given by t'char(p) = dpj. One finds smooth 
segments of the density which, depending on the initial data, may evolve into shocks. These 
are density discontinuities where the density jumps from a value pieit{xs,t') to Pright{xs,t') 
at the shock position Xs- Shocks travel with velocity Vg = (jVight — jicft)/(Pright — Pieft)- 
and are stable if the Lax condition fchar(picft) > Vs > fchar(Pright) is satisfied. An initial 
density discontinuity which does not satisfy the stability condition for shocks evolves into a 
rarefaction wave which is a smooth entropy solution of the continuity equation ([1]) {4]. 

On the microscopic scale, a shock is a sharp increase of the local density, averaged over a 
finite lattice segment. The shock position has diffusive fluctuations around its deterministic 
mean displacement Xs{t') — Xs{tQ) = Vgit' — tg) |2l|. The microscopic objects corresponding 
to the characteristics are spatially localized finite perturbations of the local density which 
travel with the collective velocity fcou = fchar of the lattice gas [22|]. These microscopic 
perturbations are analogous to kinematic waves [2_33 appearing in the macroscopic PDE- 
description of nonequilibrium many body systems. In view of this correspondence, we shall 
apply the intuitively appealing term kinematic wave (in slight abuse of language) also to 
travelling microscopic perturbations observable on the lattice scale. 
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In this way one can compute the macroscopic evolution of an initial density profile pro- 
vided that the initial density po{x) is subcritical everywhere. For hopping rates Wn which do 
not lead to condensation (pc = oo) this is not a restriction. Furthermore, the open system be- 



haves in a way analogous to usual lattice gases with open boundaries 



24, 



25 



26 
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29|. 



One can compute the evolving density profile from ([T]) even if t 
is met, since the bulk density remains subcritical at all times 



le condition for condensation 



30|. 



III. HYDRODYNAMICS IN THE CONDENSATION REGIME 

The standard considerations of the previous section fail for po{x) > pc- For po{x) > pc 
the current cannot be computed from the product measure and already at po{x) = pc the 
interpretation of the characteristics as local perturbations away from the local density be- 
comes open to doubt. A fluctuation below pc would surely travel with collective velocity 
■^^coii = Vchar, but the interpretation of a fluctuation above pc becomes dubious. In order 
to argue that nevertheless the continuity equation ([T]), properly interpreted, describes the 
macroscopic time evolution of the density under Eulerian scaling we now consider a su- 
percritical density segment po{x) > pc- The key observation is that local stationary is a 
sufficient, but not a necessary condition for ([1]) to be valid. 

To fix ideas we first assume an initial distribution of particles such that the whole chain 
has supercritical density pk > Pc- In a finite system with particles the stationary current 
does not increase beyond 0c in the thermodynamic limit L, N oo with p = N/ L fixed even 
though finite-size corrections can be substantial {s]. Therefore the current-density relation 
takes the form (Fig. [1]) 

{6 for < p < Pc 
(pc for p > Pc. 

This behaviour can be rationalized by viewing the site ko where the condensate is located 
in the infinite system as a source of particles which are emitted at constant rate Woo onto 
site ko + 1. On the other hand, site k^ acts as a sink which absorbs all particles arriving 
from site k^ — 1. The condensate site itself contains always an infinite number of particles, 
defined such that for every finite L one has N = pL particles in the chain. This effectively 
breaks the ring into an open chain with a source at the left boundary where particles are 
injected with rate Woo and a right boundary (at L oo) where particles are absorbed. In 
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FIG. 1: Schematic picture of the stationary current density relation in the ZRP with condensation. 
Above the critical density the current is constant with jc = 4>c- 

this case one has indeed j = (f)c = jc 0, Q- 

For densities not too far above pc, there is an initial metastable regime of finite dura- 
tion [11]. After this the system starts to coarsen, i.e., "small" condensates at finite distance 
start to evolve such that larger condensates grow at the expense of smaller condensates. In 
.he asyooet* ZRP .he ,„ean distance ,,o.s as flS. ™tn eventually only a single 
very large condensate remains. This condensate moves on the lattice on time scales where 
p > 1 depends on the details of the hopping rates We conclude that on time scales 
t' = Lt the mean separation between condensates is proportional to Hence the local 
density p{x,t'), coarse-grained over length segments proportional to L, remains unchanged 
on the Eulerian time scale. On the other hand, between the condensates the system has an 
average background density pc and the current has its stationary value j = jc irrespective 
of p jol. Therefore dxj{p) = and ([1]) is satisfied even though the system is not stationary, 
but coarsens. 

Furthermore, one expects that kinematic waves travel with critical collective velocity 
= fcoii(p*) between condensates and that they are absorbed into a condensate when 
they reach it after a time of order \/L. Hence the long-time average collective velocity 
vanishes, thus allowing for identification of fcou with f^har = dpj = for p > pc- 

It remains to investigate the situation where a finite fraction of the chain is initially 
subcritical, whereas some neighbouring domain, also of length proportional to L, is super- 
critical. Inside each domain the dynamics are described by the previous considerations. In 
order to investigate the boundary between the domains we consider first a domain boundary 
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characterized by pfc(O) < Pc for k < ko and pfc(O) > pc for k > ko. We consider two piecewise 
constant density profiles with pieft (pright) as the density of the subcritical (supercritical) 
domain. Without loss of generality we set kg = 0. Hence the two domains are connected by 
a density jump which on the macroscopic scale corresponds to a shock. In the supercritical 
domain A; > one expects at Eulerian time scale a series of small condensates, separated by 
a fiuctuating background with average density pc. We define kMt = 0{VL) as the position 
of the leftmost condensate. In the whole supercritical domain the current is jc- 

In the subcritical domain (pieft < Pc) one has j < jc- Hence the infiux into the su- 
percritical domain is less than the fiux inside and as a result the domain boundary moves 
towards the first condensate. When it reaches fcieft the fiux j from the subcritical region 
to its left into the condensate becomes smaller than the fiux jc out of the condensate into 
the supercritical region. As a result the condensate shrinks in size and finally disappears. 
Then the domain boundary moves on until it hits the second condensate and the process 
of condensate annihilation sets in again. Due to mass conservation the domain boundary 
thus moves into the supercritical domain and "eats it up". The shock separating the two 
domains is stable according to the usual stability criterion since ficft > Vg > "Wright- Notice 
that, as discussed above, fright = for Pnght > Pe- 
lf Pioft = Pc, mass conservation results in a vanishing shock velocity. Microscopically one 
expects a coarsening domain with background density pc to the right of the critical domain 
which has p = pc, but no condensates. The microscopic structure of the shock is not a jump 
in the background density, instead it originates from the condensates in the supercritical 
domain. The microscopic position of the shock is determined by the fluctuations of the 
position of the leftmost condensate. 

Now we consider the space-reflected case where pieft > Pc and Pright < Pc- In this case 
the rightmost condensate in the supercritical region serves as a source with constant flux 
jc that feeds into the subcritical domain with j < jc- Thus one expects a constant density 
proflle with p = pc to the right of the rightmost condensate up to the beginning of the 
subcritical domain. Hence effectively one has a critical region (initially of a size of the order 
a/L) connected to a macroscopic subcritical domain. The domain boundary moves with 
collective velocity v*^ii and to its right a rarefaction wave develops according to the entropy 
solution of ([1]). For v*^^^ = the critical domain with pc does not grow on the Eulerian time 
scale. 
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Therefore, for any domain boundary between subcritical and supercritical segments, the 
coarse grained time evolution of the density profile can be computed from ([1]) with the 
prescription 

{dpi for < p < Pc , . 

(8) 
for p > Pc. 

The characteristic velocity of the hydrodynamic equation and the collective velocity of the 
lattice gas coincide. Notice that for the TAZRP j = (p cind therefore for p < Pc the collective 
velocity is proportional to the inverse of the compressibility 

«: = 0|^ = ^((iV^)-(iV)') (9) 
of the lattice gas in the grand-canonical ensemble. 

IV. SIMULATION RESULTS 

In this section we present the results of Monte Carlo investigations which confirm the 
preceding theoretical analysis. We simulated a model in which particles hop to the right 
with rates 







Wn=l + -. (10) 

n 



In [61] it was observed that, for this choice of Wn, one sees condensation for b > 2. The 
critical density pc = 1/(6 — 2) and the critical collective velocity is given by [Oj 

for 2 < 6 < 3 

-'coll - ^ , (11) 

'-^^^^ for6>3. 

Note that the collective velocity vanishes for 2 < 6 < 3 because the compressibility is infinite. 

In Figs. [2] and [3] we show results for b = 2.5 and 6 = 4 respectively. In both cases 
we simulated a periodic lattice starting from an initial condition of Pk{0) = 2pc for < 
k < L/2 and Pfe(O) = elsewhere. The lefthand columns of the figures contain snapshots 
of the coarse-grained density profile for a single realization at increasing times whereas the 
righthand columns show an additional average over stochastic histories. At the left boundary 
of the supercritical region one clearly sees a right-moving shock front which "eats up" the 
condensation regime. Figure [3] also clearly demonstrates the growth of a critical domain 
(p = Pc) whose right boundary moves with a speed consistent with (fTTj) . 
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(a) t = 0.4 X 10'*, single realization 
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(c) t = 0.8 X lO"*, single realization 
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(e) t = 1.6 X 10"*, single realization 
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(b) t = 0.4 X 10'*, average over histories 
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(d) t = 0.8 X 10^, average over histories 
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(f) t = 1.6 X 10^, average over histories 



FIG. 2: Monte Carlo simulation results for TAZRP on a ring of size L = 2 x 10^ where Wn is given 
by (jlOp with b = 2.5. The data points are coarse-grained density profiles, averaged over 5 x 10^ 
sites, at increasing times for a single realization (left column) and an average over 10 realizations 
(right column). The initial profile (solid line) and final density (dotted) are shown for comparison. 
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The ensemble average gives qualitatively the same picture as the coarse-grained space 
average for a single realization, thus illustrating the self-averaging nature of the process. 
For longer times, diffusive fluctuations smooth out the the sharp domain boundaries when 
the profile is averaged over histories. 



V. CONCLUSIONS 



In summary, we have argued that the hydrodynamic description ([T]) of the ZRP remains 
valid for supercritical densities even though locally the system is not stationary under Eu- 
lerian scaling. In order give a meaning to ([1]) we have analyzed the coarsening process and 
found that the continuity equation has to be supplemented by the results ([7]), (jHj) for the 
current. We have demonstrated the validity of the theory by Monte Carlo simulation of the 
TAZRP. 

We expect similar analysis to be valid for the symmetric ZRP, under diffusive scaling 
L ^ oo with k = xL, t' = tL'^. In this case one gets the conservation law dt'P + d'^(j){p) = 0. 
Inside a supercritical domain one has 0(p) = (pc, corresponding to a vanishing collective 
diffusion coefficient. Adapting the considerations of the coarsening process to this situation 
one arrives at a solution that is analogous to the free boundary solution of the phase- 
segregation problem in the low-temperature phase of higher dimensional lattice gases 2l|. 
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(a) t = 0.4 X 10"*, single realization 
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(c) t = 0.8 X lO'', single realization 




■ 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

(e) t — 1.6 X lO"*, single realization 
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(b) t = 0.4 X 10'^, average over histories 
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(d) t = 0.8 X 10^, average over histories 
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(f) t — 1.6 X 10"*, average over histories 



FIG. 3: Same as Fig.[2]but for 6 = 4. From (f) we estimate the speed of movement for the righthand 
bomidary of the critical domain (i.e., the flat section with p ^ Pc = 0.5) as 7500/16000 ~ 0.47, to 
be compared with the theoretical prediction, from (llip . of 4/9 ~ 0.44. 
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